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In this paper, we propose an in-depth review of the vector and tensor polarizabilities of the two 
energy levels of the 87 Sr clock transition whose measurement was reported in [P. G. Westergaard et 
al, Phys. Rev. Lett. 106,210801 (2011)]. We conduct a theoretical calculation that reproduces the 
measured coefficients. In addition, we detail the experimental conditions used for their measurement 
in two Sr optical lattice clocks, and exhibit the quadratic behaviour of the vector and tensor shifts 
with the depth of the trapping potential and evaluate their impact on the accuracy of the clock. 

PACS numbers: 06.30.Ft, 42.62.Fi, 37.10.Jk 


Introduction 


Optical lattice clocks are now the most stable fre¬ 
quency references m, and their accuracy is steadily 
improving towards 10 -17 to 10" 18 [Ml- In these clocks 
the systematic effects due to the atomic motion are can¬ 
celled by trapping a large number of ultra-cold atoms 
(typically 10 4 ) in the Lamb-Dicke regime, using an op¬ 
tical lattice formed by a standing wave laser beam. A 
specificity of these optical clocks is the strong light shift 
induced by the intense trapping light. At the so-called 
“magic wavelength”, for which the scalar polarizabilities 
of the fundamental and excited clock states are identical, 
this light shift is largely suppressed [7J- However, in or¬ 
der to achieve a high accuracy, the residual polarization- 
dependent and higher order light shifts have to be eval¬ 
uated. In reference [5] , we reported exhaustive measure¬ 
ments of these light shifts with two 87 Sr lattice clocks, 
and showed that they are compatible with an accuracy 
of 10~ 17 at a trapping potential of 100 recoil energies. We 
reported the first experimental resolution of the vector-, 
tensor-, and hyper-polarizability and put an upper bound 
on higher order effects for 87 Sr. 

In this paper, we focus on the first-order electric dipole 
interaction to propose a theoretical study of the coeffi¬ 
cient of the decomposition of its Hamiltonian in vector 
and tensor irreducible operators that can reproduce the 
experimental results reported in [5]. The first section 
introduces the irreducible operator formalism describing 
the atomic polarizability. In the second section, we make 
use of this decomposition to theoretically evaluate the 
vector and tensor polarizability coefficients. The last 
section details the measurement of these coefficients and 
show their agreement with the theoretical estimates. Fi¬ 
nally, we describe the non-linear dependence of the vector 
and tensor shifts with the depth of the trapping poten¬ 
tial. Furthermore, we report on precise measurements of 
the magic wavelength and its sensitivity that provides in¬ 
sight in the physical properties of electronic levels of 87 Sr 
that are useful for characterizing the clock accuracy [5], 


I. FACTORIZATION OF THE 
POLARIZABILITY OPERATOR 

The Hamiltonian that describes an atomic level \(f>) = 
| nJFm) in the presence of the electromagnetic field of a 
trapping light and a static bias magnetic field B s reads: 

H = Hz + H eS . (1) 

Hz is the Hamiltonian representing the Zeeman interac¬ 
tion: 

H z = 9J ^P-B s , (2) 

n 

where g is the Lande factor and ps the Bohr magne¬ 
ton. H e ff is the effective Hamiltonian representing the 
the off-resonant electric dipole interaction between the 
atom and the trapping light of complex amplitude £ or 
power density T = Eqc \£\ 2 /2 in second-order perturba¬ 
tion theory m- 

H eS =--a |£| 2 , (3) 

where the polarizability operator a reads: 

^ = 1 £* • d\W) W\z ’ d | £■ d\4>'){^' |e* ■ d 

ft Wo — W; W 0 + U>1 

I r) 

where the sum runs over all excited states \(f>') = 

| n' J'F'm'). d, is the electric dipole operator, Huj o is the 
energy difference between | <//) and | </>), and u>i and e are 
the angular frequency and polarization of the trapping 
light, respectively. Because the two states involved in 
the clock transition (5s 2 1 So and 5s5p 3 Pq, as show in 
figure [l]) both have J = 0, their respective polarizabili¬ 
ties do not a priori depend on the magnetic sub-level m 
considered. However, the hyperfine Hamiltonian describ¬ 
ing odd Sr isotopes with a non-zero nuclear spin slightly 
breaks the rotational invariance and introduces a minute 
dependence of the polarizabilities of 1 S'o and 3 Pq on m 
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given by equation ([5|. The first term of this equation, 
for j = 0, is the scalar polarizability. It is the main, rota- 
tionally invariant contribution to the polarizability. The 
vector polarizability (j = 1) appears only when F > 1/2 
and when the light has a non-linear polarization. It is 
equivalent to a fictitious magnetic field along the light 
wave vector. It is an odd function of to, and accordingly 
vanishes when the light shift is averaged over opposite 
values of m. The last term (j = 2), appearing if F > 1, 
is the tensor polarizability that results in a polarization 
dependent polarizability with an even dependence in to. 

The following section is devoted to the theoretical cal¬ 
culation of the three polarizability coefficients for the two 
atomic levels involved in the 87 Sr clock transition. 


FIG. 1: Energy levels for Sr. The two clocks states are rep¬ 
resented in red. The F = 9/2 hyperhne component of the 
three green states are involved in the vector and tensor po¬ 
larizabilities of the excited clock state 5s5p 3 Po, as explained 
in section El The other states are considered to calculate the 
polarizability coefficients of the clock transition. 


and t. This dependence can be made analytically ex¬ 
plicit by expanding the polarizability operator as the sum 
of three irreducible operators m EH- I n the (2 F + 1)- 
dimensional basis |m) = {| nJFm), —F < m < F} of the 
hyperhne manifold of the Hilbert space describing either 
clock state, the matrix elements of this expansion read: 


II. THEORETICAL CALCULATION OF THE 
POLARIZABILITIES 


In order to further calculate the polarizabilities, the kj 
coefficients can be related to the \<j>') —> \<j>) transition 
rate T^>, or to the radiative lifetime r of the excited 
state \<j>') through: 


KF'lldHF)! 2 _ Q 3ire 0 Hc 3 r _ D 3n£ 0 Hc 3 u ( 

-^hfs q dxb' -^hfs q -^fs •> 


2F' + 1 


Wn 




where the fine and hyperhne branching ratios read: 


(9) 


<TO 2 |dK) = ^^ 

3 =0 M=-j 

(5) 

where the first sum runs over the tensor rank j. In this 
expression, the dependence in £ is exclusively contained 
in the tensor product [t (g> t*]jM, the dependence in to 
in the Clebsch-Gordan coefficient _ M and the de¬ 

pendence in | (j)') in the coefficients: 


I n'J'F'\ 


1 


(-l) j 


u>0 — UJi LOq + LOi 


+ 


( 6 ) 


where kj is expressed using Wigner 6j-symbols and re¬ 
duced dipole elements: 

(7) 

The canonical scalar a s , vector a v and tensor a t polar¬ 
izabilities are then dehned by rescaling these coefficients: 


cx s \/3 
and at = 


2 F 


F + l 

I 2F(2F - 1) 

3(F + l)(2F + 3) 


K i, 


K,. 


( 8 ) 


Bis = (2J +1){2L' + 1)^ l s J L , | , (10) 
Bhis = (2F+1)(2J , + 1)|^' \ J,j , (11) 

and £ is a dimensionless parameter accounting for the 
hne structure splitting. 

When neglecting the hyperhne interaction, the two 
states of the clock transition with J = 0 do not exhibit 
vector or tensor polarizabilities, as can be seen by setting 
F = 0 in equation (1 t|. However, the hyperhne splitting 
of the excited states |p'): 


Gifs 


\ AC B IC(C + !)-/(/ + 1)J'{J' + 1) 
2 2/(2/ — 1)J'{2J' — 1) 


( 12 ) 


where C = F'{F' + 1) - /(/ + 1) - J'{J' + 1) and A and 
B are magnetic dipole and electric quadrupole constants, 
gives rise to vector and tensor components. The relative 
magnitude of these latter components is expected to be 
on the order of the ratio between the hyperhne splitting 
i^hf s and the frequency of the optical transitions |</>) —» 

W) E3j: 


&v,t ^ ^hfs 
&s ^optical 


(13) 


These three physical parameters are sufficient to com- However, for a state \<f>) with J = 0, we can show that for 
pletely explicit the matrix of the polarizability operator any value of /, and with J' = 1 according to the electric 
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dipole selection rules: 


i+J' 

Y_] k 2 C = 0 


(J = 0 , J' = 1), (14) 


F'=I-J' 


such that the contribution of the magnetic dipole term 
A to the tensor polarizability is cancelled at first order 
in j/hfs- Hence, the dominant contribution to the tensor 
shift is the quadrupole B , or the second order in v\ x f s . 
whichever is the largest. No such cancellation occurs for 
the vector shift. This property is particularly observed 
for 8 ' Sr, as the excited states with the largest hyperfine 
splitting feature like 5s6s 3 Si and 5s5p 3 Pi. As 

a consequence, their contribution to the tensor polariz¬ 
abilities of the clock states is orders of magnitude smaller 
than what could be expected from (131. 


Using experimental measurements of the lifetimes of 
excited states, completed by theoretical estimates of the 
strength of atomic transitions Mi we can calculate the 
dynamic polarizabilities for various states of 87 Sr at the 
magic wavelength (w; = 27r x 368.6 THz). 

The scalar, vector and tensor polarizabilities of the fun¬ 
damental state 5s 2 1 tS'o are derived by summing over the 
excited states 5s5p 3 P± and 5skp 1 Pi with k = 5... 20. 
They read, expressed in atomic units (a.u.): 


5s 2 1 S’o 

a s 

OLy 

ot t 

279.8 

4.75 x 10“ 5 

1.57 x 10 -5 


The main contribution to the scalar polarizability comes 
from the 5s5p X P\ state, while the vector and tensor terms 
are dominated by the 5s5p 1 Pi and 5s5p 3 Pi states. As 
expected, the ratio between the vector (resp. tensor) po¬ 
larizability and the scalar polarizability is on the order of 
10“', comparable with the ratio between magnetic dipole 
term A - about 200 MHz (resp. the quadrupole term B - 
about 50 MHz) and optical frequencies - about 500 THz. 

The polarizability of the 5s5p 3 Pi(P = 9/2) is esti¬ 
mated by summing over the bsks 3 S\ states with k = 
6... 10, the 5 p 23 P states and the 5 skd 3 D states with 
fc = 4... 9: 


5s5p 3 Pi(F = 9/2) 

Us 

OLy 

u t 

304.4 

-50.3 

-101.2 


As expected from the fact that J ^ 0 for 3 P \, the vec¬ 
tor and tensor polarizabilities are on the same order of 
magnitude as the scalar polarizability. Figure[2]confronts 
these numerical values to the experimental spectroscopy 
of the 5s 2 1 5 0 —¥ 5s5p 3 Pi(F = 9/2) inter-combination 
line. 

In the case of the excited clock state 5s5p 3 Po, a di¬ 
rect calculation is not sufficient because state mixing 
with the 5s5p 3 Pi(P = 9/2), 5s5p 3 P 2 {F = 9/2) and 
5s5p 1 Pi(F = 9/2) states alters its physical properties. 



Detuning 6v (MHz) 


FIG. 2: Spectroscopy of the saturated 5s 21 S’o —> 
5s5p 3 Pi(P = 9/2) transition by depletion of 87 Sr atoms in 
an optical lattice, for a linearly polarized trapping light, no 
magnetic field, and a trap depth of 1520 Er, where Er = 
h x 3.47 kHz is the recoil energy associated with the absorp¬ 
tion of a lattice photon with frequency wi. In this configura¬ 
tion, the polarizability operator |5| is already diagonal when 
choosing the quantization axis along the lattice polarization. 
Its eigen-values are a(m) = a a + ctt 3 "jq~2F < Lr)~ 1) ■ The graph is 
centred on the light shift free frequency of the transition, and 
displays the splitting between magnetic sub-levels induced by 
the scalar and tensor polarizabilities of both states. The ar¬ 
rows indicate the theoretical light shift 8v = — l«( m )\£\ 2 for 
each value of \m\ , as calculated from equations (15 I and (161. 


Since these states have 4 / 0, their vector and tensor po¬ 
larizabilities are large and may contribute to the polar¬ 
izabilities of the 3 P 0 state. This state mixing is written: 

10) ='52 c p\ ( I , p)o, ( 17 ) 

V 


where the coefficients c p are US]: 


\4 > p)o 

5s5p 3 P 0 

5s5p 3 Pi 
(P = 9/2) 

5s5p 3 P 2 
(P = 9/2) 

5s5p 1 P 1 
(P = 9/2) 

Cp 

1. 

2.3 x 10" 4 

-1.38 x 10 -6 

-4.1 x 10“ 6 


(18) 


When state mixing is involved, expanding the matrix 
elements of the polarizability operator d4j) is more in¬ 
volving, but the result is that equations ([5]) through © 
remain valid expressions of these matrix elements, pro¬ 
vided they are averaged over all (|</> p )o, |</>q)o) combina¬ 
tions with a weight c p c*, and provided the branching ra¬ 
tios are replaced by: 


B ls = (_i)Sp+s 9 +2£'+A-K/ ? ^/(2j p + i)(2 J q + 1) 


x 


( 22 / + 1 ) 


J' 1 J p 
LSI’ 


J' 1 J q \ 
L S L' J ’ 


(19) 


P\\ fs 


(2U+1)(2J'+1) 


F' 1 F 
J p I J' 


F' 1 f) 
J q I J'}' 


( 20 ) 
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Using these expressions, we can show that for J p = 0 and 
J q = 1, ^2=0 such that this combination 

has a negligible contribution to the tensor shift, and for 
J p = 0 and J q = 2, Y^F'=\i-j'\ = 0 such that this 

second combination has a negligible contribution to the 
vector shift. 

Among all the (|</> p )o, |<A?)o) combinations involving the 
states listed in the table of equation (18), only a few lead 
to a non-negligible contribution. They are listed in bold 
characters in the following table: 


5s5p 3 P 0 

I^p)o 

<Aj)o 

a s 

OL-u 

a t 

3 P 0 

3 P 0 

288.8 

9.3 x 10- 3 

1.9 x KT 6 

3 P1 

3 Pi 

2 x 10" 5 

2.7 x 10" 6 

-5.4 x 10 6 

3 Po 

3 Pi 

2 x 10" 6 

1.82 x lO’ 1 

-1.5 x 10~ 6 

3 Po 

3 p 2 

< 10" 13 

-1.8 x 10" 8 

3.70 x 10 4 

Total 

288.8 

1.91 x 10" 1 

3.65 x 10" 4 


( 21 ) 

In this table, the (|0 P ) O , |0 9 )o) and (\4> q )o, \4> P )o) contri¬ 
butions for p ^ q are grouped together. As expected from 
the theoretical considerations above, the tensor term of 
the ( 3 P 0 , 3 Po) configuration is negligible, and the main 
contribution to the vector and tensor polarizabilities re¬ 
spectively come from the mixing combinations ( 3 Po, 3 Pi) 
and ( 3 P 0 , 3 P 2 ). 

To express these polarizabilities in units more adapted 
to experimental purposes, we note Uq = \a s \S | 2 the trap 
depth. At the magic wavelength, a s and hence Uq are by 
definition identical for the two clock states. The numer¬ 
ical values of a s shown in equations (15) and (21) show 
that this fact is approximately rendered by our calcu¬ 
lation, which also gives an order of magnitude of a few 
percent for its accuracy, arising from uncertainties on the 
decay rates from excited states. For the vector and tensor 
polarizabilities, we introduce the coefficients: 


ft 


V 


a v 1 1 
a s h 2 F 


and ft* 


a t 1 1 

h 2F(2F - 1)' 


( 22 ) 


The following table gathers the theoretical values for 
k v and K* for the two clock states, as well as their differ¬ 
ence Aft u,t = ft"’*( 3 P 0 ) - «"’*(%): 


10) 

ft 1 ’ (mHz/F/j) 

ft* (pHz/E R ) 

3 P 0 

-255 

-60.9 

% 

-6.5 x 10" 2 

-2.7 

3 Po - 'So 

-255 

-58.2 


III. EXPERIMENTAL DETERMINATION OF 
THE POLARIZABILITIES 

During the clock operation, a static bias magnetic field 
is applied to split the magnetic sub-levels. The energy 




FIG. 3: Left: parametrization of the light polarization e and 
the static bias magnetic field B s . e x (resp. e y ) is a unit 
vector along the major (resp. minor) axis of the trapping 
light polarization e, such that the polarization ellipsis (red) 
is included in the (e x ,e y ) plane. Right: projection in the 
plane defined by the trapping light wave vector k and B s . 
The fictitious magnetic field representing the vector light shift 
adds to B a to form an effective magnetic Held b e s. If the 
vector component of the Hamiltonian dominates its tensor 
component, the Hamiltonian is diagonal for a quantization 
axis ez defined along 6 e ff- When the polarization e*is elliptical, 
the quantization axis therefore depends on the trap depth Uo, 
introducing non-linear effects in the Zeeman shift and average 
clock frequency as a function of Uo- 


levels are then the eigenvalues of the full Hamiltonian 

written in equation (llj). Because F is a vector operator, 
the Zeeman Hamiltonian adds up to the vector compo¬ 
nent of the polarizability operator. Therefore, the vector 
component of the Hamiltonian is equivalent to the action 
of an effective magnetic 5 e ff field equal to the vector sum 
of the static magnetic field B s and the fictitious magnetic 
field that represents the vector light shift (see figure [3]) : 


UfF 


K v i{ef\ £*)Uo , 5 

I s 

9 


(24) 


where g = gp-s/h. If the quantization axis ez is chosen 
along the effective magnetic field (i.e ez = & e ff/||^eff||), 
the vector component of the Hamiltonian is diagonal 
in the basis {| m), — F < m < F} of eigen-vectors of Fz- 
However, the tensor part of the Hamiltonian is a priori 
not diagonal in this basis. Yet, since k v ft*, if we as¬ 
sume that the bias field is large enough for the Zeeman 
splitting to be much larger than the energy splitting due 
to the tensor shift (For a trap depth Uo ~ 100 Er, this 
assumption is satisfied if B s 1 pT), the states {|to)} 
remain approximate eigen-states of the full Hamiltonian. 
Expanding the tensor product in equation (|5j) then yields 
the energy shift 5v of level | m) of a given hyperfine man¬ 
ifold: 


5v(m) =+mg\\b e s\\+iJpUo, (25) 


with 


P = (3|e ■ e z \ 2 ^ l)(3m 2 - F(F + 1)). (26) 

The middle term of equation ( [25] ) includes both the 
vector light shift and the Zeeman shift. It can be isolated 
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from the scalar and tensor terms by measuring the half 
difference Z s {m) = \ {8v{rn) — 5v(—m)) between oppo¬ 
site magnetic sub-levels: 

Z s = ra<?||& e ff|| = m\J ( gB s cos 9 + k v £U 0 ) 2 + ( gB s sin df, 

(27) 

where £ = 11 i(eA £*) 11 is the degree of circular polarization 
of the trapping light, and 9 is the angle between the wave 
vector k and the bias field B s (as shown on figure[3]). The 
Taylor expansion of the later equation to second order in 
Uq yields: 

( k v £) 2 

Z s ~ mgB s + mn v t; cos 9 Uq + m sin 2 9-— —— Uq. (28) 



To experimentally evaluate the differential vector polar¬ 
izability coefficient Ak v , we measured A Z s = Z s (^Pq) — 
Z.fSo) for different trapping depths with m = ±9/2 and 
a circular polarization for the trapping light (|£| = 1). 
From these data we extrapolated the derivative of the 
Zeeman shift at zero trap depth dAZ s /dUo(Uo = 0) = 
mAn v £ cos 9. The trap depth is measured by observ¬ 
ing the longitudinal motional side-bands of the trapped 
atoms m- their spacing and shape give a direct mea¬ 
surement of the trap depth (at the trap centre) and of lon¬ 
gitudinal and transverse temperatures of the atoms PH- 
From these three parameters, we can deduce the average 
trap depth experienced by the atoms. 

Repeating the experiment for various values of 9 en¬ 
ables to deduce the differential coefficient, as reported 

in 0: 


\Ak v \ = 0.22 ± 0.05 Uz/E r , (29) 


in agreement with the theoretical estimate reported in 
equation (23). The uncertainty is limited by the knowl¬ 
edge of the polarization state e in the lattice cavity. In 
the usual clock operation, B s is orthogonal to the wave 
vector k ( i.e. cos 9 — 0). This configuration minimizes 


the first order term of equation (28), but it also maxi¬ 


mizes its quadratic term. This latter term can therefore 
easily be observed, even for moderate trapping depths 
below 200 Er. Figure [4] shows such a quadratic de¬ 
pendence for an elliptical polarization. In order to fit 
the data points of this figure, the two Lande factors of 
3 P 0 and 1 5'o have to be known. We performed a precise 
measurement of the ratio between these coefficients by 
locking the clock on and n transitions, as shown on 
figure [5} 


g( 3 Po) 

g^So) 


= 1.58794(7). 


(30) 


Given the Lande factor of the fundamental state g( 1 S'o) = 
184.4(1) [151 US] , this gives: 

g( 3 P 0 ) = 2.928(2)MHz/T, and A g = 1.0842(7)MHz/T. 

(31) 


FIG. 4: Differential Zeeman shift A Z s of the 87 Sr clock res¬ 
onance as a function of the trap depth Uo for an elliptical 
polarization |£| ~ 0.9 and |m = 9 /2|. The experimental 
data points are fit with equation (271, assuming k v ( 1 So) = 0. 
These four data points are sufficient to determine the three 
free parameters of this equation that describe the geome¬ 
try of the experimental setup, here B s = 355.0 ± 0.2 /rT, 
£k v ( 3 Po) = 212 ± 2 mHz /Er, and 9 = 63.9 ± 0.5 mrad. 



Bias field (|_iT) 


FIG. 5: Ratio R between the differential Zeeman shifts A Z a 
(for linear polarization, so that no dependence in Uq is ob¬ 
served) for cr(|m| = 9/2 —> \m\ = 7/2) and 7r(|m| = 9/2 —> 
\m\ = 9/2) transitions. This ratio is related to the ratio be¬ 
tween Lande factors by g( 3 Pq)/^ 1 Sq) = 9 (R — 1)/(9 R — 7). 


The measurement of the tensor shift was conducted 
with a purely linear polarization for the trapping light, 
and therefore with a quantization axis along the bias field 
B s . The tensor component is extracted by measuring the 
average clock frequency shift for various values of /?. The 
results are reported in table [T] 

a _ v( 3p o, m) - v(}S 0 ,m) + v( 3 P 0 , -m) - p^Sq, -to) 

LW — - 

2 

= (Ak s + PAk^Uq. (32) 


in agreement with previously reported values m- 


where: 
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Meas. # 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

sin# 

1 

1 

0.92 

0.92 

1 

1 

1 

1 

0.99 

0.99 

0.955 

0.955 

1 

COS 0 

0.602 

0.602 

0.602 

0.602 

0 

0 

1 

1 

1 

1 

1 

1 

0 

m 

9/2 

7/2 

9/2 

7/2 

9/2 

7/2 

9/2 

7/2 

9/2 

7/2 

9/2 

7/2 

9/2 

p 

3.12 

1.04 

-2.96 

-0.99 

-36 

-12 

72 

24 

70.0 

23.3 

63.05 

21.0 

-36 


TABLE I: Geometrical configurations chosen to measure the tensor shift. A strictly linear polarization is selected by the 
enhancement cavity of the optical lattice by using an intra-cavity birefringent element. For this polarization, the /? coefficient 
reduces to j3 = (3 sin 2 6 cos 2 0 — l)(3m 2 —T'(T’-l-l)) = (3 cos 2 —l)a(3m 2 — F(F+1)) where a is the angle between the polarization 
and the bias field B s . We switch between cos 0=1 and cos 0 = 0 by changing the polarization eigen-mode of the lattice cavity. 
The angles are derived from the relative values of the Zeeman shift when the magnetic field is changed. 


Aa s _ a s ( 3 Pp) - a a ( 1 S 0 ) 

Qt s 


(33) 


is the remaining differential scalar light shift due to a 
possible detuning of the trapping light from the magic 
wavelength. The resulting value for the tensor shift co¬ 
efficient is then, as reported in j§j: 

A A 

A«* xp = = (-57.7 ± 2.3) /xHz /E r , (34) 


in agreement with the theoretical value shown in equa¬ 
tion ( |23[ ) well within the experimental error bar. 

We now consider the behaviour of the tensor light shift 
when the light polarization e is not perfectly linear. In 
this case, the quantization axis e% for the excited state 
3 P 0 is aligned with the effective magnetic field b e # and 
thus varies with the trap depth Uq. Consequently, the /3 
coefficient depends on Uq through: 


£-ez\~ = 


£■ B s 


IMI 


- \£-£b s 


, K v ^ COS 0 Uq /tt ^ 
1-2 gB. + ° m 


(35) 

where ep s is the unit vector along B s . Because of this 
coupling effect between the vector and tensor light shifts, 
the average clock frequency contains a frequency shift 
A u vt quadratic with the trapping potential reading: 


A v vt = —"f v t\£- eB„| 2 £cosf?[/g, 


with 


7«t = (3m 2 - F(F + 1)) 


6K t ( 3 P 0 )K v ( 3 P 0 ) 
~g( 3 P 0 )B s ’ 


(36) 

(37) 


assuming that exhibits no vector shift. For 87 Sr we 
have 7 vt = 11 fjMz/E^ for B s =0.1 mT. This quadratic 
frequency shift on the average clock frequency is three 
orders of magnitude smaller than the quadratic term of 
the Zeeman shift written in equation ( |28[ ). But this 
effect is quite significant when compared to the hyper¬ 
polarizability coefficient 7 = 0.4 /rHz[SITS]. How¬ 
ever, it is dramatically reduced when the polarization is 
linear, B s _L k and B s _L e. Its experimental demonstra¬ 
tion is for now challenging, as it is predominant in a con¬ 
figuration where the vector light shift is large and blurs 
the atomic resonances. 


Finally, repeated measurements of the total light shift 
over two years with our two strontium clocks [5] con¬ 
firmed the experimental value for the magic wavelength 
published in [ 8 j, for which the differential scalar light shift 
is cancelled (i.e. Ak s = 0): 

u magic = 368 554 725 ± 3 MHz. (38) 


The uncertainty on this value is limited by our knowledge 
of (3 used to subtract the tensor light shift from the total 
light shift from equation (32). In addition, we measured 
the sensitivity of the scalar light shift with the trapping 
light frequency Z 7 att : 


3Ak s 


-15.5 ± 1.1 ^Hz/Pr/MHz. 


(39) 


This value is in agreement with the estimation given in [9j 
from Monte-Carlo simulation using transition strengths 
and scalar polarizabilities of the clock states at various 
frequencies. 


Conclusion 

In this paper, we have conducted a theoretical esti¬ 
mate of the vector and tensor polarizabilities of the 87 Sr 
clock states in agreement with experimental measure¬ 
ments. We have shown that the differential vector po¬ 
larizability is largely due to state mixing of the excited 
clock state 3 P 0 with 3 P\. The tensor polarizability is 
mainly due to state mixing of 3 Po with 3 P 2 , with a small 
contribution from the hyperfine structure of 1 P\ on ^o- 
We also described a non-linear behaviour of these light 
shifts that are largely cancelled in normal clock opera¬ 
tion, but that have to be considered as the accuracy of 
optical lattice clocks continues to improve. 
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